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ABSTRACT 



■ Aims. The circumstellar disk of AB Aurigae has garnered strong attention owing to the apparent existence of spirals at a relatively 
young stage and also the asymmetric disk traced in thermal dust emission. However, the physical conditions of the spirals are still not 
well understood. The origin of the asymmetric thermal emission is unclear. 

^ . Methods. We observed the disk at 230 GHz (1.3 mm) in both the continuum and the spectral line ^^CO J=2^1 with IRAM 30-m, the 

O Plateau de Bure interferometer, and the Submillimeter Array to sample all spatial scales from 0'.'37 to about 50''. To combine the data 
obtained from these telescopes, several methods and calibration issues were checked and discussed. 
^ , Results. The 1.3 mm continuum (dust) emission is resolved into inner disk and outer ring. The emission from the dust ring is highly 

O I asymmetric in azimuth, with intensity variations by a factor 3. Molecular gas at high velocities traced by the CO line is detected next 

I ■ to the stellar location. The inclination angle of the disk is found to decrease toward the center. On a larger scale, based on the intensity 

O ' weighted dispersion and the integrated intensity map of ^^CO J=2^1, four spirals are identified, where two of them are also detected 

I in the near infrared. The total gas mass of the 4 spirals (Mspirai) is 10~^ < Mspirai < 10~^ M©, which is 3 orders of magnitude smaller 

C/^ . than the mass of the gas ring. Surprisingly, the CO gas inside the spiral is apparently counter-rotating with respect to the CO disk, and 

■ it only exhibits small radial motion. 
Conclusions. The wide gap, the warped disk, and the asymmetric dust ring suggest that there is an undetected companion with a 
mass of 0.03 M© at a radius of 45 AU. The diff'erent spirals would, however, require multiple perturbing bodies. While viable from an 
energetic point of view, this mechanism cannot explain the apparent counter-rotation of the gas in the spirals. Although an hypothetical 
fly-by cannot be ruled out, the most likely explanation of the AB Aurigae system may be inhomogeneous accretion well above or 
below the main disk plane from the remnant envelope, which can explain both the rotation and large-scale motions detected with the 
30-m image. 
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1 . Introduction later evolutionary stages. As a result, the structure and kinemat- 

^. tt t. 1 . n 1 1 . ics of the disks of pre-main- sequence stars provide clues to the 

• • Circumstellar disks are comnaon y found around pre -main- ^.^j^ ^^^^^^^ ^^^^ ^^^^^ formation. 

. ^ sequence stars (IStrom et al.lll989l: iBeckwith et al]ll990h . They ^ ^ 

K> are believed to play an important role in processing the ex- With the improvement in angular resolution and in sensitiv- 

^ cess angular momentum of cloud cores and in allowing the ity, CO maps reveal that these disks usually exhibit Keplerian 

C3 ; accretion to proceed inward toward the central star radially kinematics on scales of 100s AU (e.g. Pietu et al. 2007). In 

via g ravitational instabilities (see the review bv iDurisen et alJ the meantime, more and more circumstellar disks are found to 

'2007) or magnetorotational instability (Balbus & Hawlev"l99l'; have complex structures, such as the large cavities foundin 

Balbus & Hawley 1998), or vertically (Konigl & Salmeron HD 135344B (Brown et al. 2009), MWC 758 (Isella et al. 20T^), 

201 1, and references therein). Within these disks , planetary sys- LkCa 15 (Pietu et al. 2006), and HH 30 (Guilloteau et al. 2008). 

tems are expected to form in the early (Inutsuka et al.ll201Q>) or A more thorough analysis of the structure and kinematics of 

circumstellar disks are therefore important for providing con- 

* Based on observations carried out with the IRAM 30-m telescope straints on the physical conditions of planetary system forma- 

and Plateau de Bure interferometer. IRAM is supported by INSU/CNRS tion. 
(France), MPG (Germany), and IGN (Spain). 

- Based on observations carried out with the Submillimeter Array located at about 140 pc, AB Aurigae (hereafter AB Aur) 

(SMA). The SMA is a joint project between the Smithsonian one of the closer Herbig Ae stars of spectral type of 

Astrophysical Observatory and the Academia Sinica Institute of AO-Al (Hernandez et al. 2004). Mod eling of its mid-mf rared 

Astronomy and Astrophysics and is funded by the Smithsonian spectral energy distrib ution (SED) by iMeeus et alJ (l200ll) and 

Institution and the Academia Sinica. iBouwman et al.l (|200Q|) shows that AB Aur belongs to Group I, 
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where the SED can be reconstructed by the combination of a 
power law and a black body, so this young star is usually con- 
sidered as the prototype of the Herbig Ae star surrounded by a 
large flaring disk. 

Indeed, AB Aur is surrounded by a large amount of mate- 
rial that is o bserved from the n ear infrared (NIR; e.g. HST-STIS 
imaging by Gradv et al."l999') up to the mm domain ( Pietu et al. 
^05; Corder et al. 2005). The NIR images have revealed a flat- 
tened reflection nebula, clos e to pole-on, which ex tends up to 
r ~ 1300 AU from the star. Ipuka gawa et al.l (l2004l) also made 
deep NIR images of AB Aur using the Coronographic Imager 
Adaptive Optics systems on the Subaru telescope. Their obser- 
vations revealed spiral features within the circumstellar matter. 
Interestingly, part of these spirals ha ve also been mar ginally de- 
tected in CO 3-2 line emission by iLin et al" using the 
Submillimeter Array (SMA). Finally, independent studies of 
the CO line kinematics have reveal ed that the outer d isk struc- 
ture i s in sub-Keplerian rotation (iPietu et al.l l2005l: iLin et al.l 
[20Q6h. unlike what ha s been found in other similar systems 
(e.g. lSimon etan '2000'). Since the disk does not appear massive 
enough to be self-gravitating, one plausible explanation would 
be the youth of the syste m. The Platea u de Bure interferometer 
(PdBI) images from Pietu et al. l (l2005l) at 1.3mm and in CO 2-1 
also show that the material is truncated toward the center with a 
ring-like distribution with an inner radius around 70 - 100 AU. 
The origin of such a cavity has not yet been identified. This may 
be due to a substellar or a very low-mass stellar companion lo- 
cated at ~ 20 - 40 AU from the primary, as suggested by the 
recent very large array (VLA) obse rvations at the wavelengths 
of 3.6 cm by Rodriguez et aL] (l2007h . 

Therefore, both the gas and dust orbiting around AB 
Aur present several puzzling features, which include spirals, 
(sub)Keplerian motions, and an inner cavity. This object appears 
very diff'erent from the simple picture of a flaring, Keplerian gas 
and dust disk. It may represent the real complexity of an evolv- 
ing young system in which planet formation has already started. 
In this paper we attempt to improve our physical understanding 
of AB Aur by conducting a joint IRAM-SMA project in order to 
better define the mm/submm properties of the source. Sections 2 
and 3 present the observations and results. We discuss the impli- 
cations of the observed properties in Sect. 4 and summarize in 
Sect. 5. 



2. Observations and data analysis 

We imaged the 1.3 mm continuum emission by utilizing the 
highest angular resolution (0'.'3) data from the PdBI, together 
with the short-spacing information from the SMA. The gas 
kinematics were studied with ^^CO J=2^1 emission using the 
IRAM 30-m, PdBI, and SMA. 

2.1. IRAM 30-m data 

The IRAM 30-m ^^CO J=2^1 data were obtained on Sep 9- 
10, 2009. The receivers were tuned in single sideband, with a 
measured sideband rejection of about 13 dB. We used a spectral 
resolution of 40 kHz (0.052 km s"^) for the spectra. We used the 
on- the- fly mapping technique. A fixed position near AB Aur was 
used as the reference position to subtract the background. This 
position was, however, not free of CO emission: it was observed 
in frequency switching to recover the full flux, and the spectrum 
was added back to the data obtained in the on-the-fly processing. 
Consequently, the 30-m is not missing any spectral emission. 



However, the continuum level is lost due to the linear baseline, 
which was fitted to remove residual atmospheric contamination. 

Pointing calibration was done in two steps. First, pointing 
was checked every two hours at the telescope. Focus was also 
monitored and found to be stable. Pointing accuracy was then 
refined by obtaining a series of independent images, each taking 
about half an hour of elapsed time. Each image was then recen- 
tered (through cross-correlation) on the nominal position deter- 
mined from the PdBI data, using the high-velocity channels that 
only come from the disk emission (see later discussion). This al- 
lows us to correct the 30-m pointing errors as a function of time. 
Finally, each spectrum was corrected for the determined pointing 
off'sets (using linear interpolation in time from the values derived 
for each map), and the final images were produced. 

The overall calibration procedure guarantees a registration 
better than V between the 30-m and interferometer data, and an 
amplitude calibration precision better than 10 %. 

2.2. Plateau de Bure interferometer mosaic 

The high angular resolution data were obtained in mosaic cen- 
tered on AB Aur, observed with long baselines in the AB con- 
figuration, on Jan 30, Feb 20, and Mar 14, 2009, and on Jan 22 
and Feb 2, 2010, in both continuum and in ^^CO J=2^1. The 
seven-field mosaic uses a hexagonal pattern with field separated 
by IV\ about half of the half-power primary beam of the 15-m 
dishes at this frequency (23''). The baselines range from 45 m to 
760 m. These newly obtained data were suppleme nted by the ex - 
isting compact configuration (CD) data from Pietu et al.l (l2005l) . 
which cover one field (central field) and baseline lengths from 
15 to 174 m. 

In all cases, the correlator was set with two units covering the 
^^CO J=2^1 line at high spectral resolution (80 kHz, or 0.104 
km s"^). The remaining six units were used to observe the con- 
tinuum emission covering a bandwidth of 0.96 GHz in HH and 
VV polarization. All data were processed using the CLIC soft- 
ware in the GILDAS package. Phase and amplitude calibration 
were performed using the quasars 3C111 and 0518-1-134, while 
the flux calibration used MWC 349 as a reference. 

The combined data set results in an angular resolution of 
a.'52 X a'37 at positional angle (PA) of 24° and a.'56 x a.'42 
at PA of 27° for the 1.3 mm continuum and ^^CO J=2^1, re- 
spectively. The noise level of the 1.3 mm continuum, (Xi, is 0.23 
mJy beam"^. The noise level for the ^^CO J=2^1 line is 22 mJy 
beam"^ (2.2 K) at the spectral resolution of 0.203 km s"\ which 
is the resolution of the maps shown later in this paper. 

2.3. SIVIA observations and data reduction 

The PdBI mosaic lacks short spacings below 45 m except for the 
central field, and the IRAM 30-m can only provide reliable spac- 
ings up to 15 m. Furthermore, the 30-m data do not give any con- 
tinuum information. We therefore observed in ^^CO J=2^ 1 and 
in continuum at 1.3 mm, a single field with the SMA in the sub- 
compact configuration, to sample the deprojected baseline range 
6 to 68 m (4.4 to 52 kA). The correlator was set to cover the 
^^CO J=2^1 line in the upper sideband with high spectral reso- 
lution (50.8 kHz, or 0.066 km s"^). The rest of the chunks with 
the equivalent bandwidth of 3.28 GHz were set for the contin- 
uum. 

The short spacing data from the SMA were obtained on Sep 
14, 2009 under weather conditions of r = 0.1 at 225 GHz and 
with humidity between 7% and 20%. The data were calibrated 
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Table 1. Observation parameters 



Parameters 



Obs. scheme 



RMS 
(mJy/beam) (K) 



Angular resolution 
sizerx") PAC) 



30-m CO 



on-the-fly 



1150 



0.2 



12x12 



PdBI 



1.3 mm 
CO 



seven-field 
mosaic 



0.23 
22 



0.027 
2.2 



0.52x0.37 
0.56x0.42 



24 
27 



SMA 



1.3 mm 
CO 



single 
pointing 



1 

110 



8x10" 
0.09 



7.4x4.4 
7.4x4.4 



166 
166 



and reduced using the MIR software. MIR is a software pack- 
age to reduce data taken with the SMA. The bandpass, flux, and 
gain calibrator are 3C454.3, Uranus, and 3C111, respectively. 
The phase scatter on 3C111 is ± 20°. The angular resolution of 
the SMA dataset is 7! A x 474 at PA of 166°, and the noise level is 
1 mJy beam"^. The half-power field of view of the SMA at this 
frequency is 55'', which is sufi&cient to cover the field covered 
with the PdBI with a seven-field mosaic. 

2.4. Data merging 

All spectral data were resampled to a 0.203 km s"^ spectral res- 
olution. Corrections for the proper motion of AB Aur, (1.7, - 
24.0) mas yr"^ from Hipparcos, were applied to all data, and all 
positions referred to equinox J2000. From the existing observa- 
tions, several combined datasets were produced. In this paper, 
we present and discuss the following cases. 

1. The combined 30-m + SMA uv data: they were imaged in 
a standard way. The resulting data cube has a typical rms 
noise of 0.13 Jy beam"^ in nearly empty channels. In the 
channels with strong emission, the rms rises up to 0.39 Jy 
beam"\ as a result of the eflTective dynamic range limitations 
due to phase (SMA data) or pointing (30-m) errors, as well 
as uncertainties in intensity calibrations. The peak signal-to- 
noise ratio is about 160. 

2. The PdBI mosaic: seven-field mosaic observed with 
medium, long, and very long baselines (up to 760 m) but 
lacking baselines shorter than 45 m. The central field has a 
muc h better uv coverag e because of the existing PdBI data 
from lPietuetaP (l20Q5h . 

3. The PdB+SMA+30-m data: the two previous data sets were 
combined using a specially designed procedure. To do so, we 
utilized the fact that the SMA+30-m deconvolved image has 
high signal-to-noise. We thus used it as a model for the sky 
emission to compute model visibilities and added these to the 
PdBI uv table for each field. In this process, the respective 
primary beams of the SMA (centered) and PdBI (diflTerent 
for each field) must be taken into account. In the combina- 
tion, the choice of the uv coverage on which visibilities must 
be computed is rather arbitrary. At the minimum, we can use 
the sampling defined by the SMA tracks completed with the 
30-m spacings up to about 10 m. For the spacings beyond 
10 m the SMA becomes more precise than the 30-m owing 
to the limited pointing accuracy of the latter. Alternatively, 
we can use any arbitrary sampling up to about 40 m, beyond 
which the SMA+30-m is no longer sensitive enough com- 
pared to the PdBI data. The weights must also be adjusted to 
match the PdBI weights and to obtain a final distribution of 
weights that is as uniform as possible; this reweighting nec- 
essarily implies some loss in sensitivity. This overall process 
is strictly equivalent to considering that the SMA+30-m de- 
convolved image (corrected for the primary beam) is like a 
30-m image, except for the diflferent angular resolution. We 



experimented with several weightings and uv samplings, and 
as expected, the results are relatively insensitive to the (arbi- 
trary) choices, provided we do not select uv spacings that are 
too extreme from each telescope. 
4. In case the signal-to-noise ratio in the SMA+30-m image is 
not high enough, a diflTerent procedure may be used. Rather 
than retaining the complete image, we can instead use the 
Clean components, as only noise should be left out. Each 
Clean component must be corrected for the diflTerent pri- 
mary beam attenuations (including the pointing center oflT- 
sets for the PdBI). A uv table can then be created from these 
Clean components. We have, however, enough sensitivity in 
all cases that this procedure gives results that are totally com- 
patible with the one mentioned in item (3) above. 

Relative calibration of the data from these three telescopes 
was checked by comparing visibilities from the CO channels (the 
continuum flux being too weak for this purpose). This process 
did not reveal any major discrepancy, but had limited accuracy 
(20 % level at best) because of the complex spatial structures. 
Comparison of strictly redundant spacings does not provide suf- 
ficient sensitivity. We thus kept the original flux calibration scale 
of all data sets, which, given the observing techniques, should be 
more accurate than 10 %. In practice, our conclusions are fairly 
insensitive to the relative calibration errors, since the only com- 
bined data being used is for the continuum (see Sect l3.1.l1 to 

The maps and spectra of ^^CO J=2^1 line presented in this 
paper are from PdBI data, unless it is explicitly mentioned. All 
relative positions cited in this paper are with respect to the coor- 
dinate {a, ^)j2ooo = (04^55'^45r84, +30°33'04r30). Table 1 sum- 
marizes the observational parameters. 

3. Results and modeling 

3.1. Continuum observations 

Figure [T] is a montage of our 1.3 mm continuum maps obtained 
from the SMA and PdBI data. 

3.1.1. Large-scale continuum 

For the SMA data (upper-left panel in Fig. [T]), the total flux is 
about 110 mJy and the rms noise 1 mJy beam"^ This value is 
only marginally lower than the single-dis h measurement of 13 
mJy, made at slightly higher frequency ("Henning et al.l ll998h . 
Although the emission appears localized toward AB Aur, it does 
not only originate in the circumstellar disk of AB Aur. The PdBI 
data only recovers about 80 mJy of the total flux (see below). 
After removing the ring emission recovered by the PdBI from 
the SMA data, there is residual emission toward the south of the 
disk (upper-right panel of Fig.[T]). 

The excess emission peaks about V!4 southeast of AB Aur, 
and half of the missing flux is within one SMA beam. There is 
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_SMA continuum 




_ SMA continuum with 
subtraction 





Offset (arcsec) 




J 

-20 




-PdBI continuum with 
. -. ring subtraction 




30°33'06" ^ 



30°33'05" - 




Offset (arcsec) 



Fig. 1. Continuum emission maps at 1 .3 mm of the AB Aur region. Top panels: Larger scale image obtained from the SMA data (left 
panel) and the residual image after removing the emission detected by the PdBI (right panel). Contour levels are -2, 2, 4, 6, 8, 10 
to 80 by 10 X 1 mJy per 7V4x4V4 beam. Bottom panels: Higher angular resolution images obtained with the PdBI. Left panel: The 
blue, red, and black pluses mark the peaks of CO 2-1 highest blue-shifted, red-shifted peak, and the stellar location, respectively. 
The black cross marks the northeast peak of the 1.3 mm continuum. The dotted contours mark the intensity of 2 (Xi, where (Ti= 0.23 
mJy per 0'.'52xO'.'37 beam. Black contours start from 4 cri in steps of 2 (Ti. Middle panel: Residual emission from the PdBI after 
removal of a best-fit smooth ring. The inner and outer radii are indicated. Contour levels are multiples of 2 (Ti. Right panel: the solid 
white ellipse marks the 1.3 mm peak ring with a radius of 145 AU, / of 23° and PA of -30° centered on the 1.3 mm stellar peak. The 
blue and red contours mark the CO 2-1 emission starting from and stepping in 66 mJy per 0'.'56x0'.'42 beam at the velocity channels 
of 2 and 9.3 km s"^ respectively. The 3.6 cm continuum emission are in thick black contours at the level 3, 6, 9, 12, 15 x 9.3 yuJy 
per O'.' 27x0^.^26 beam. The synthesized beam of 1.3 mm and 3.6 cm are labeled as black and white ellipses in the lower left corner, 
respectively. 



no significant structure on scales smaller than about 3'\ Such 
structures would have been recovered by the PdBI, because it 
provides suflftcient uv coverage and sensitivity to 2 - 3'' scales. 
This excess emission is thus smooth on a scale of 3'' or more. It 
presumably traces the remnant envelope of AB Aur. 

With a dust absorption coeflftcient, k = 0.02 cm^ g"^ 
gas-Fdust) at 230 GHz, which is characteristic of disks 
Beckwith et al.l ll990), and a dust temperature of 30 K, the to- 
tal gas mass obtained with the SMA would be 10"^ Mq. For 
more pristine, ISM-like dust, the mass would be about four times 
higher. For the residual emission, the mean A^h is 1.3 x 10^^ 
(cm"^). Assuming the residual emission is spherically symmet- 
ric, the mean number volume density is 5 x 10^ (cm"^). 

3.1 .2. High angular resolution data 

At high angular resolution, the 1.3 mm continuum emission is 
resolved into three structures (bottom panels in Fig.[T]): 

- a continuum peak at the central stellar position, 

- a continuum gap between the central peak and a dust ring. 



- a dust ring/ellipse. 

The continuum peak near the center is detected at the level 
of 5.7 (Ti with a total flux of 1.3 mJy. It is not resolved within the 
0^.^37 (50 AU) interferometric beam. The mm continuum peak 
also coincides (within the astrometric accuracy of our measure- 
ments, about 0^.^05) with the stellar position and the 3.6 cm emis- 
sion peak (Fig.[T] lower-right p anel). The flux density of the 3.6 
cm emission is 0.2 mJy (Rodri guez et al.L 2007). With a spectral 
index of 0.6, characteristic o f optically thick emission in an ex- 
panding ionized jet or wind (iPanagia & Felli|[l975) , this would 
yield 1.5 mJy at 1.3 mm, so that the 1.3 mm signal from the star 
position could be purely due to an ionized jet. On the other hand, 
optically thin free-free emission would have an essentially flat 
spectrum, so the remaining 1 . 1 mJy flux could come from dust 
emission from an inner disk. A very small (2-3 AU), optically 
thick, and dust disk would provide enough flux. A (gas-hdust) 
mass of 4xlO"^Mo is required for such inner disk, assuming the 
temperature of 80 K and k = 0.02 cm^ g"^ This mass is likely a 
lower limit since the dust emission should be partially optically 
thick. 
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The continuum gap (cavity) between this central peak and 
the dust ring has a width of ~(y/6, corresponding to about 90 
AU. It appears more pronounced in the east and southeast, but 
this is not statistically significant. 

The dust ri ng/ellipse confirms the structure found by 
iPietu et aP (I2QQ5 ). With higher spatial resolution, we are able 
to further separate the emission from the central star. For an 
assumed dust temperature T(r) = SOCr/lOOAU)^ "^ K (see the 
following paragraph) and k = 0.02 cm^ g"\ the ring mass is 
5.0 ± 0.5 X 10"^Mo. With the same assumptions as in Sect. 3.1.1, 
the peak A^h and nu of the dust ring are 7.3x10^^ (cm"^) and 
7.0x10^^ (cm"^), respectively. The intensity contrast on the dust 
ring is < 3. 

The continuum emission is further fitted using DiskFit 
with an inclined, sha rp-edged ring, as done for GGTau 
(iGuilloteau et al.|[r999h . We fitted two diff'erent functionals: a 
truncated power law given by 

^(r) = ^o(r/Ror' for Rin < r < R^uu (1) 

and an exponentially tapered distribution given by 

^(r) = ^o(r/Ror'^M-(r/Rcf-n for R,n < r. (2) 

The best-fit of the dust ring parameters are listed in Table[2l Both 
functionals give similar results, in particular for the inner radius. 
The derived surface density profiles are given in Fig. [21 The ta- 
pered edge solution deserves specific comment. With the large 
inner radius also found in our best-fit solution, the high negative 
value found for the exponent is required to fit the steep decrease 
in surface brightness in the outer part. A solution in which both 
the inner and outer parts are fit by a profile with no sharp inner 
radius is possible, but would imply even more negative expo- 
nent values , below -4.5, a situation very similar to that found by 
llsella et al.l (l2010h for MWC 758. Besides being somewhat bet- 
ter than this solution, our fit demonstrates that the outer edge is 
rather sharp and that the dust emission does not extend beyond 
about 250 AU. 



Table 2. Fitting parameters of dust ring 





Power law 


Exponential tapering 


Inner radius (AU) 


108 + 6 


109 + 4 


Outer radius (AU) 


193+4 


[350] 


Critical radius (AU) 


[oo] 


157 + 4 


P 


0.8 + 0.5 


-2.4+1.0 


PA(°) 


-31 + 6 


-27 + 9 


in 


20 + 3 


22 + 4 



13.3b . The intensity of this northeast peak is 0.76 mJy beam"^ (3.3 
(Ti). Because it is statistically not significant, we do not discuss 
this component further. 

The center of the dust ring is at (a,^)J2ooo=(4^55'^45.85^ 
30°33'4V32), consistent with the stellar location. The position 
angle (PA) of the rotation axis with respect to north of the best- 
fit ring is — 30°. The inclination angle, /, is -23° based on the 
ratio of minor axis to major axis assuming the ellipse is thin. 
After subtracting the ring, the residual image given in Fig. [T] 
(bottom-middle panel) shows a strong northsouth residual on the 
western side and a strong negative residual on the eastern side, 
suggesting that the emission on the dust ring is not azimuthally 
symmetric. 

3.2. CO gas 

The channel maps and the spectrum obtained with the IRAM 
30-m toward AB Aur are shown in the top panel in Fig. [3l The 
angular resolution of the IRAM 30-m is 1 V!2. The emission is 
extended and dominated by the envelope near the systemic ve- 
locity. There is a velocity gradient in the southeast-northwest di- 
rection. The combined image of the 30-m -h SMA CO data with 
the smoothed resolution of 7f4 is presented in the lower panel 
in Fig. [3] The peak flux density is 84 Jy beam"\ which corre- 
sponds to 40 K. At this resolution, the emission toward the south 
is clearly seen near the velocity of 6.9 to 7.3 km s"^ . 




50 100 500 

Radius (AU) 

Fig. 2. Best-fit surface density distributions for the power-law 
(blue) and tapered-edge (red) models. Dotted curves represent 
the ± Icr of the best-fit value. Without an inner edge, the tapered 
edge solution would provide too much flux from the inner 100 
AU, as shown by the dashed line. 



In addition to these three structures, there is a northeast con- 
tinuum peak 2! 2 away from the stellar positi on, corresponding 
to the ET peak identified bv lLin et al.1 (l2QQ6h in their 0.88 mm 
continuum image. It is in between CO SI and S2 arms (see Sect. 



3.2.1 . Large-scale map at high angular resolution 

We present in Fig.|4]the ^^CO J=2^ 1 mosaic image from PdBI. 
The full data set was also studied, combining the 30-m data with 
the SMA field and the PdBI mosaic. A comparison between the 
30-m-hSMA-hPdB recombined image and the PdB-only mosaic 
reveals that all the interesting detected structures are actually 
confined to within the central beam of the PdBI mosaic, with 
few significant structures (at high resolution) beyond about 12'' 
from AB Aur. Moreover, the extended emission recovered in the 
30-m-i-SMA-i-PdB image makes the structures less obvious. The 
PdB-only mosaic channel maps presented in Fig. [4] are more 
suitable to pinpointing the contrasts between diff'erent structures. 
Since these images do not include the extended flux, the appar- 
ent brightness temperatures are lower limits, and correction for 
the missing flux must, however, be included to derive physical 
parameters, such as temperature. 



3.2.2. Inner gas 

The ^^CO J=2^1 emission is detected from the velocity of 1.6 
km s"^ to 9.5 km s"^ (Figs. |4]and[5]). The emission peaks in 
the high- velocity ends (Fig. [5]) are near the 1.3 mm continuum 
central peak. The mean position of the highest velocity CO peak 
is ofl'by ~ 0^.^05 to the 1.3 mm stellar peak. Relative positions are 
controlled by the bandpass calibration accuracy, which would 
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Fig. 3. Top panel: ^^CO J=2^1 channel maps from the 30-m observations at ir/2 resolution. The field of view is 90'\ The insert 
shows the spectrum toward AB Aur at higher spectral resolution. Contour spacing is 2 K. Bottom panel: combined ^^CO 30-m + 
SMA image at 7.4'' resolution. Contour levels are 1 to 5 K by 1 K, then every 5 K (1 K = 0.42 Jy beam"^ at this resolution). In the 
cleaning, the residual are around 0.2 Jy beam"^ in the 16 inner channels. 
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Fig. 4. Channel maps of the ^^CO J=2^1 line from the PdBI mosaic. The contours are 1 to 9 by 1 x 0.1 Jy per Of! 56 x Of! 42 beam 
(4.5(T, 9.91 K). The Vlsr is marked in each panel in km s"^ . The red cross marks the stellar location. Cyan arcs mark the CO spirals 
identified in Sect. 3.3.1. The S4 arm is a dotted arc, because it does not have as broad a linewidth as SI, S2, and S3. 



■1.21 • ■ 

; y: ; 


1.62 


2.02 


2.43 


2.84^^^^ ■ 


.' ' u 
■ 8.7 


5 


■ 9.14 


■ 9.54 


9.95 


■10.36 ■ 

: ■% ] 


• 













04^55"'45'9 45'8 



R.A. 

Fig. 5. The same as Fig. 3 but toward high- velocity channels. 
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Table 3. Best-fit parameters of high velocity CO emission 



Fixed / (°) 


[23] 


[29] 


Inner radius (AU) 


< 20 


< 20 


T (100 AU) (K) 


54 + 4 


54 + 4 


q 


0.20 + 0.08 


0.14 + 0.07 


PA(°) 


-38.2+1.1 


-38.2 + 1.1 


Vrot(lOOAU) (km.s-1) 


5.60 + 0.15 


4.57 + 0.12 


V 


0.50 + 0.04 


0.72 + 0.03 


M, (Mo) 


3.5 + 0.15 


[2.35 + 0.12]* 



Notes. Best-fit parameters are derived assuming two different inclina- 
tion angles: i=23° from best fit to dust emission, and i=29° from stellar 
properties and Keplerian assumptions. Vrot is the deprojected rotation 
velocity at 100 AU. Values in brackets are fixed. 
Stellar mass only meaningful if v = 0.5. 



yield relative positions valid to about 1/100* of the beamsize, 
and by the signal-to-noise. The diff'erence is at 2 to 3 cr level and 
is only marginally significant. 

High- velocity CO gas (Fig.O is detected well within the dust 
ring, at least as close as 20 AU from the star. This high- velocity 
emission is consistent with a rotating disk with the rotation axis 
of PA = -38° ± 1°. However, because of the restricted velocity 
range (more than 2.4 km s"^ away from the systemic velocity) 
and limited spatial extent, the inclination angle is poorly con- 
strained. 

Physical quantities based on two inclination angles are fur- 
ther derived and presented in Tabled We use d a simple disk 
model and the DiskFit code (Pietueta D 12007') to derive them: 
the temperature, T, and velocity, V, are assumed to follow a 
power-law dependence on the radius, r: T = To(r/ro)~^ and 
V = Vo(r/ro)~^. On one hand, if the inclination angle is fixed to 



23°, corresponding to the best fit of the continuum ring, the best- 
fit velocity law will be compatible with Keplerian rotation. The 
stellar mass (M^) will be 3.5 ± 0.2 Mq, which is higher than the 
mass of 2.4 Mq expected from the AB Aur spectral type. On the 
other hand, if the stellar mass is fixed to be 2.4 Mq and the disk 
is assumed to be in Keplerian rotation, the best-fit inclination an- 
gle will be 29°, which deviates from the best-fit value for the dust 
ring of 23°. Finally, if the inclination angle is fixed to be 29°, the 
best-fit velocity power-law index becomes v = 0.72 ± 0.03, i.e., 
the Keplerian rotation assumption is rejected at the Icr level. We 
do not find a common model that can fit both the dust ring and 
high- velocity CO emission. 

In comparison, on a larger scale, the rotation axis as derived 
from ^^CO by Pietu et al. (2005) has the PA = -30° ±1°, which 
is consistent with the best-fit value of the dust ring. The incli- 
nation angle is in between 36° and 42°, which is significantly 
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Fig. 6. Maps of m oments 0, 1, and 2 from left to right panel for ^^CO J=2^1 (upper row; this work) and ^^CO J=2^1 (lower row; 
iPietu et al.ll2QQ5l) . The white ellipse and arcs trace the best-fit dust ring and the identified CO spiral arms, respectively. Left panels: 
The contours start from and step in 0.12 Jy beam"^ km s"^. Middle panels: Contours steps are 0.4 km s"^ Right panels: Contours 
steps are 0.1, 0.2, 0.3, 0.4, 0.6, 0.8, 1.0, ... km s"^. The unit of the wedge is Jy beam"^ km s"^ for moment maps and is km s"^ for 
moment 1 and moment 2 maps. 



larger than the values derived from our new high angular res- 
olution data. We note that the velocity field on a large scale is 
sub-Keplerian. Another puzzling issue is the derived systemic 
velocity, Vsys. It is 5.73 ± 0.02 km s"^ as derived from the high- 
velocity parts only and is 5.85 ± 0.02 km s"^ for the outer disk. 
This diff'erence is marginally significant at the level of 4cr. Such 
a deviation in Vsys may be due to an asymmetry in brightness 
between the blue- and red-shifted wings, reflecting that the dis- 
tribution of the gas is inhomogeneous. 



3.2.3. The disk 

The integrated intensity (moment 0), the intensity weighted 
velocity (moment 1), and the intensity weighted linewidth 
(moment 2) ma ps of ^^CQ J=2 ^1 from this work and of 
^^CO J=2^ 1 by Pietu et al.1 (l2005h are shown in Fig.[6l As a ref- 
erence, a white solid ellipse that marks the 1.3 mm dust ring is 
also indicated. Figure [6] shows that ^^CO emission is much more 
extended (7'' or 1000 AU) than ^^CO, which is strongly domi- 
nated by the rotating disk, and extends out to 4'' (560 AU). The 
outer radius in continuum is even smaller, at most 280 AU (see 
Fig.O. While opacity difl'erences naturally produce this ordering 
(R(dust) < R(^^CO) < R(i2C0)), the sharpness of the dust profile 
would lead to only marginally larger radii for ^^CO and ^^CO if 
opacities only were to explain this result. O ur result thus con- 
firms and reinforces the proposition made by lPietu et al.l (l2005h 



of a change in dust properties (emissivity or gas-to-dust ratio) 
at a radius around 280 AU. Our more precise determination is 
now in e xcellen t agreement with the feature detected at 20 yum 
by Pa ntin et al.l (|2005). This suggests that a change in dust prop- 
erties is more likely than a drop in gas-to-dust ratio beyond 280 
AU. 

3.3. Spiral structure 

3.3.1 . Identification of the spiral arms 

The ^^CO moment m aps exhibit very d iff^erent morphology from 
those of i^CO (from |Pietuetal.ll20Q5k . While ^^CO traces the 
rotating disk, there are spiral-like patches with broader veloc- 
ity dispersion and excess integrated intensity compared to the 
Keplerian patterns in the ^^CO line. 

The ^^CO moment and moment 2 maps from Fig. [6] allow 
us to identify four "spiral" patterns, CO SI to CO S4. From the 
moment 2 image, there are three patches where the linewidths 
are broader. We marked them as CO SI to CO S3. Part of 
these arms can also be seen in the moment map. Although 
the linewidth of CO S4 is not as broad as CO SI to CO S3, there 
is excess emission at S4 in the moment map. Besides the 0'.'5 
maps, we also checked the ^^CO J=2^1 maps with different 
tapering. The excess emission near the spirals can still be seen 
but is not as obvious as in the 0:'5 maps. On larger scales, the 
^^CO line might be too contaminated by the emission from the 
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Fig. 7. Locations on the spirals where the spectra in Fig. [8] are 
taken. The 1.3 mm dust ring is marked as a black ellipse. 



envelope in the combined maps. Another reason might be that 
the excess emission comes from structures that are actually very 
compact and not as bright as the Keplerian disk emission. We 



note that the labeled spirals are traced by first marking the lo- 
cations where the moment and moment 2 maps clearly show 
excess emission "by eye". Then, we fit four parabolic arcs to the 
marked locations. 

We further checked the spiral emissions in the channel maps 
(Fig. [4]). The SI and S2 arms are on the blue-shifted part of the 
disk. In SI, the spiral emission can be seen at 3.9 to 6.7 km 
s"^ The S2 arm is the most extended spiral among these four 
arms: it appears from 4.7 to 7.3 km s"^ and connects to the SE 
part of the disk. The S3 and S4 arms are on the red- shifted part 
and at the near side of the disk. The S3 arm has the broadest 
linewidth and it shows up from 4.3 to 7.3 km s"^ The S4 arm 
appears from 5.7 to 7.7 km s"^ . The bases of the four spiral arms 
are diff'erent. The SI arm seems to trace back to the 1.3 mm 
gap (cavity; see the channels 3.9 to 5.5 km s"^ in Fig. [4] and 
moment 2 in Fig. (6]). The S3 arm can be traced back to the ring 
inner edge (see channels 4.3 to 4.9 km s"^ in Fig.|4]). S2 and S4 
arms originated somewhere near the dust ring. The connection of 
these CO spirals with the dust ring and the NIR spirals and the 
comparison with the spirals detected in CO 3-2 will be discussed 
in Sect ion 

Spectra of ^^CO J=2^1 and ^^CO J=2^1 line along the 
spirals at locations spaced every Off 56 (see Fig. [7] for point 
labeling) are given in Fig. [S] While the ^^CO J=2^1 line 
shows only the emission from the (nearly) Keplerian disk, the 
^^CO J=2^1 spectra clearly exhibit a new kinematic compo- 
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nent. The ^^CO J=2^1 spectra can be decomposed into two 
Gaussians, one of which traces the rotating disk (at almost the 
same velocity as the ^^CO J=2^1 line) and the other traces 
an excess emission on the spiral. They are named ^^COdisk and 
^^COspirai, respectively. Results of the Gaussian best-fit are listed 
in Table m and the kinematics are discussed in Sect. 4.2.3. The 
velocity off'sets of the two Gaussians are typically 1 to 2 km 
s"^. This provides a natural explanation of the apparent broad 
linewidths in the spiral locations: there is a bulk of gas moving 
at diff'erent velocities from the rotating disk but not because of 
turbulent motions on the spirals. We note that although the spi- 
rals (arcs) are labeled by eye and not based on quantitative crite- 
ria, the parameters derived in this way are expected to reflect the 
bulk properties of these non-Keplerian component. 



3.3.2. Temperature and opacity 

The peak temperature in the CO maps, around 80 K, and the 
average brightness temperature in the inner IT', 40 K, indicate 
temperatures ranging from 30 K outside to 100 K in the inner 
regions of the disk, confirmed by the more sophisticated disk 
analysis presented in Tabled 

The CO line opacity can be estimated from the ^^CO/^^CO 
ratio map (Fig.O. On large scales, in the velocity channels 5.22 
to 6.72 km s"^ in the southern part, and 4.97 to 5.47 km s"^ in 
the northeast of the disk, the ratio ^^CO/^^CO is over 1 (Fig.©. 
This must be due to self- absorption in ^^CO, implying optically 
thick in ^^CO with a temperature decreasing as a function of 
distance from AB Aur. 

On the identified CO spirals, the ^^CO/^^CO ratio is low (in 
fact there is no robust ^^CO detection), so the ^^CO opacity can- 
not be very high. A better constraint comes from the apparent 
brightness temperature of the emission, which ranges from 10 K 
to 30 K. Assuming the kinetic temperature in the spiral is com- 
parable to the one in the disk, this moderate brightness indicates 
opacities of 0.3 to 1, since the spirals are resolved. Higher tem- 
peratures (hence lower opacities) cannot be excluded, but this 
would require a multiline study. 



3.3.3. Gas mass and energetics inside the spirals 

Under LTE conditions, the column density of the CO spirals can 
be estimated using the standard equation 



2J-h 1 



/hBJ(J-h 1) 



(3) 



where A^coj is the column density of the CO molecule in the 
quantum state J (J+1 — » J). Here, Ncoj is given by 



^co,j = 2.07 X 10- 



; 2J+ 1 

27 + 3a3 



Jedv, 



(4) 



where B, y, Aui, and Tb dv are the rotational constant, the fre- 
quency in GHz, the spontaneous transition rate, and the bright- 
ness temperature Tb in K integrated over the linewidth dv in km 
s"\ respectively. Here, Z is the partition function, 2jlo (2J+1) 
exp[-hBJ(J -\- l)/kTex], and Tex is the excitation temperature. 

We assume Tex= 50 K as representative of the average 
temperature inside the arms. The derived column density of 
^^CO and ^^CO in the spirals are thus given by Nuqq = 
8.9x10^^ / Tb dv and Nuco = 6.07x10^^ / Tb dv, respectively. 
With the integrated line flux calculated using Tb and Av given 
in Table m the derived values of Nuqq range from 1.5x10^^ to 



1.5x10^^ cm"^. Using a CO/H2 ratio of 10"^, the total gas mass 
in these four spirals is only 2.5x10"^ M©. Opacity efl'ects could 
increase these values. 

An upper limit of the total mass of the spirals can be esti- 
mated from the ^^CO J=2^1 line, which is not detected in the 
spirals. With the 3cr noise level of 3 K per channel (0.25 km 
s"^), NuQQ is less than 10^^ cm"^, assuming the same linewidth 
as^fCO J=2^1. Assuming the ^^C/^^C ratio 77 in the local ISM 
(IWilsonlll994D . the upper limit of Nuqq and the gas mass in the 
spirals are ~10^^ cm"^ and 1.5x10"^ M©, respective ly. In com- 
parison, ^2^0 in the disk at 100 AU estimated in iPietu et aP 
(2006) is 1.4x10^^ cm"^. In the spirals, the column density is at 
least 10^ times lower than the one in the disk. 

The kinetic energy contained in the spirals in the rotating 
disk frame are calculated using the equation 



£:kin = 0.5MgasV^fl., 



(5) 



where Voff is the difl'erence between the Vlsr of the 
^^COspirai component and the Vlsr of the ^^CO component listed 
in TablelH The resulting £^kin on the spiral is 6x10^^ erg. The up- 
per limit of £'kin given from ^^CO J=2^1 is 3x10^^ erg, assum- 
ing the same Voff as the ^^COspirai component with the estimated 
upper limit of the spiral mass. 

If the gas on the spirals mainly moves on the disk plane, the 
reported energy will be a factor of six higher, assuming that the 
inclination angle is 23°. In comparison, the binding energy of 
a companion with a mass of 0.03 M© at 45 AU, is 2.4x10"^^ 
erg, which is at least four orders of magnitude greater than the 
kinetic energy contained in the spirals. A small perturbation of 
the orbit of such a companion is sufficient to release enough en- 
ergy to sustain such spirals. If emission is at a larger h and is 
optically thick, the estimated Ekin will be the lower limit. The 
angular momentum is estimated to be 10"^^ in cgs following the 
same method. 

4. Discussion 

At high angular resolution both in the NIR and the submm do- 
mains, the surroundings of AB Aur are complex. Figure [TOl is 
a montage of the 1.3 mm continuum image with CO spirals, 
the optical image from HST, and the NIR image from HiCIAO 
of Subaru teles cope. Superposed are th e spirals as determined 
from the NIR bv [Hashimoto et al.1 (l201 ll) and the CO spirals (this 
work). The HST image from lOrady et al] ('1999') on large scales 
also reveal the sp iral arms, which were later reported closer to 
the central star by lFukagawa et al.l (l2004l) from Subaru deep im- 
ages. 

In dust emission, the large inner cavity and well-defined 
dust ring appear much like a scaled-down version of the 
tid al cav ity -hring observed in the binary system GG Tau/A 
(iGuilloteauetal.ll 19991) . More recently, several SMA, CARMA, 
and IRAM PdBI observations have revealed the existence 
of such large cavities inside disks, such as HD 135344B, 
MWC758, LkCa 15, or even HH30. The material around AB 
Aur also diff'ers from similar objects in some important aspects: 

- In the large-scale CO emission, four spiral-like arms are de- 
tected. Part of these CO spirals are also seen in NIR even 
though the existence of a remaining envelope of dust, which 
aff'ects NIR image mostly, and gas (^^CO) partially contami- 
nates any detailed analysis of the material close to the star. 

- The disk is not in Keplerian rotation, and the dust ring 
is clearly asymmetric both in continuum and in the 
i^CO and ^^CO lines. 
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Fig. 9. The line ratio of ^^CO J=2^1 to ^^CO J=2^1 as a function of position in different velocities. There is some self- absorption 
(red color) in the southern parts. The ^^CO emission otherwise has < 1 optical depth. 



Table 4. Parameters of the Gaussian best-fit toward ^^CO and ^^CO 
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0.36 


7.02+0.02 


44.2 


0.55 


5.08+0.03 


23.4 


0.53 


3C 


-0.77 


-1.25 


-23 


6.97+0.03 


21.2 


0.30 


6.73+0.03 


28.2 


0.61 


5.25+0.04 


16.4 


0.28 


4A 


-1.24 


0.36 


54 


6.65+0.04 


18.9 


0.46 


6.79+0.04 


[40.0] 


0.40 


5.86+0.27 


21.0 


0.47 


4B 


-1.85 


0.47 


39 


6.47+0.02 


20.0 


0.28 


6.69+0.02 


42.9 


0.48 


5.86+0.03 


11.6 


0.14 


4C 


-2.44 


0.51 


29 


6.40+0.03 


17.3 


0.22 


6.71+0.02 


43.4 


0.42 


5.77+0.02 


19.9 


0.16 


4D 


-3.05 


0.51 


22 


6.27+0.04 


13.9 


0.17 


6.57+0.02 


32.3 


0.42 


5.62+0.03 


13.6 


0.16 


4E 


-3.64 


0.46 


17 


6.28+0.05 


9.6 


0.17 


6.41+0.03 


18.8 


0.43 


5.49+0.06 


18.6 


0.12 



Notes. Best-fit results of the spectra shown in Fig.[8l Coordinates of the spectra are listed in offset with respect to (R.A., Dec.)j2ooo = (04^55"*45'.'84, 
+30°33'04?30). Vlsr, Tb and Av is the observed velocity with respect to the local standard of rest, the brightness temperature and the line-width, 
respectively, ^pa is the position angle of the offset location with respect to the major axis of the disk (-121.3° from north), assuming an inclination 
of 23°. Values in brackets are fixed. Note that there is no robust ^^CO detection at locations CO S2E to CO S2J. 
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- The inner cavity is not fully devoid of gas. This gas is likely 
not uniform within the central hole (see Sect l3.2.2b . 

- Fitting the large-scale disk traced in ^^CO and C^^O (see 
Sect. 14.2.1]) and the dust ring separately leads to some appar- 
ent inconsistencies. In particular the inclination angles differ 
from 42° to 23° (while the typical errors are 1° to 3°). 

In this section, we discuss these points and see how they can 
be interpreted and possibly reconciled in a single picture. We 
focus first on a comparison of the mm/submm data with other 
similarly resolved observations and on the ring properties. We 
then discuss an overall scheme for the spiral features and present 
the dynamics of the AB Aur system. 

4.1. Properties of the ring and inner material 

4.1 .1 . Dust gap and inner material 

The 1.3 mm continuum peak at the stellar location might sug- 
gest the existence of, at least, one inner disk. The lower limit of 
the total gas-hdust mass at the stellar location is estimated to be 
4x10"^ M© (Sect. 13.1. "21) . Although we cannot rule out the possi- 
bility that the 1 .3 mm stellar peak comes from the optically thick 
free-free emission, an inner disk has been detected and resolved 
in the near infrared (NI R) and mid-in frared (MI D-IR) interfer- 
ometric observations by lEisner et a D(^004) and Id i Folco et al.l 
( 2009). The inner disk is also required to explain the observed 
accretion rate of 10"^ to 3xlO"^MrT. vr"| (Brittain et al.. 2007: 
iTelleschi et alll2007l: IPonehew & Brittainll201lh . 

Besides the continuum emission at the stellar location, the 
detection of "high" velocity ^^CO J=2^1 emission inside the 
dust cavity also indicates that this area is not devoid of gas. The 
gas mass estimated from the CO 2-1 line with velocity more than 
2.5 km s"^ away from Vsys is -10"^ M©. These gases are 50 AU 
away from the stellar location with an orbital period of 200 yr. 
From the orbital timescale argument, which the mass accretion 
rate toward the center cannot exceed the orbiting mass divided 
by the orbiting period, the upper limit of the accretion rate is in 
the order of 10"^ M© yr"\ consistent with the values reported in 
previous paragraph. 

4.1 .2. Ring properties 

Comparing the 1.3 mm continuum with NIR emission is not 
straightforward, because the dust disk has very high optical 
depth at NIR. The dust emission at these two wavelengths orig- 
inates in diff'erent dust layers. The NIR emission traces the 
scattered- stellar light by dust grains in the upper layers of the 
flared disk, while the 1.3 mm thermal emission, which is es- 
sentially optically thin, characterizes the bulk of the disk mass. 
Although the very inner parts of the NIR maps are occulted (e.g 
HiCIAO data: the area of the mask is seen in white in Fig. [TOl) , 
the superposition of the mm and infrared images reveals coher- 
ent behavior (Fig. 10). We find that the NIR emission peaks be- 
fore the mm ring at radius around 70 AU, then the mm emission 
peaks at radius 150 AU followed by a decrease in the NIR bright- 
ness, indicating that there are fewer stellar photons impinging on 
the disk behind the bulk of the dust emission. This behavior is 
fully consistent with pred ictions from stellar irrad iated-dust disk 
models (see for example JPuUemond et al.ll2007b . 

The AB Aur disk is not azimuthally homogeneous. The best- 
fit models derived from the 1.3 mm and the ^^CO data in Tables 
[2] and [3] assume uniform density distribution with azimuth. The 
residual image of the 1.3 mm continuum (i.e. after subtracting 



the best-fit ring; lower middle panel of Fig.[T]) clearly shows the 
inhomogeneity in both radius and azimuth. The azimuthal asym- 
metry can be quantified by Fig. [TT] where we trace the dust inten- 
sity at 1 .3 mm versus the position angle with respect to the major 
axis of the disk, ^pa (measured in the disk plane, by assuming an 
inclination of 23 °), inside the ring (peak location, Fig.[T]). The 
major axis of the disk is assumed to be -121.3° from north, since 
the disk ro t ation axis is at -31.3° from the best fit of ^^CO in 
Pietu et aD (l2QQ5h . There are two dips in the ring, one at ^pa of 
130° and the other at -150° (Fig. [TT1). The emission peak is at 
^PA of about 30°. The maximum intensity contrast on the ring is 
about 3. 

4.2. CO disc 

4.2.1 . Inclination angle 

The determined inclination angle / of AB Aur disk varies with 
scale. For the inner disk (radius > 20 AU), / is 23° to 29° based 
on the high velocity ^^CO J=2^1 line (Sect. 3.2.2; Table 2). 
With the best-fit dust ring at 1.3 mm (radius of 145 AU), / is 
22° to 24° (Sect. 3.1.2; Table 1), consistent with the inclination 
derived from the high- velocity ^^CO line. We note that the de- 
rived inclination angle from the ^^CO J=2^1 data is based on 
the high- velocity components (2.4 km s"^ off'set from the V^ys), 
which are less likely to be contaminate d by the envelope. On 
a larger scale (radius 70 AU to 600 AU), iPietu et al.l (|2005) re- 
ports / to be 36° based on the C^^O J=1^0 line. Because the 
C^^O J=1^0 line is optically thin and also less likely to be 
contaminated by the large-scale envelope, the discrepancy be- 
tween the values derived from ^^CO and C^^O suggests a phys- 
ical warp. We note that the low / toward a smaller scale has 
already been repor ted. Based on the NIR interferometric data, 
lEisner et al.1 (|2004 find that / is < 20° for an inner disk with a 
radius less than 1 AU. The / inferred from our new 1.3 mm con- 
tinuum and CO data is more consistent with the value derived 
from the NIR. 

Warp ed disks have been repor ted in the literature, e.g.. Beta 
Pictoris ( Golimowski et al.l l2006l and references therein). The 
AB Aur system is somewhat similar to the GG Tau A system. 
[Cruilloteau & Dut rey ( 2001) show the high angular resolution 
^^CO J=2^1 observations of CO gas inside the cavity with 
slightly diff'erent inclination angle from the one derived from the 
dust ring. 

4.2.2. Rotation 

The rotation velocity also varies with scale. In'Pietu et aD (l2Q05h 
and Lin et al. (2006), the CO emission was found to be in non- 
Keplerian rotation. This is partly due to the confusion with the 
remaining envelope since the ^^CO J=2^11ine and especially 
the ^^CO J=3^2 line are mostly optically thick, tracing both the 
envelope and the disk. However, obs erved with ^^ CO and C^^O, 
where confusion is negligible, Pietu et al. I (l2005h measured the 
velocity power-law exponents v of 0.42, 0.37, 0.47, and 0.82 for 
^^CO J=2^1, i^CO J=1^0, C^^O J=2^1, and ^^CO J=2^1, 
respectively . These exponents were fitted with all the detected 
emission in lPietu et al.l (l2005b . which is probably dominated by 
the larger scale disk. The new ^^CO J=2^1 data at higher an- 
gular resolution add an additional puzzle, since the velocities 
within 100 AU appear to exceed the orbital velocity around a 
~2.4 M© mass star. One might explain such high velocities with 
an increasing inclination angle at smaller radius, but the im- 
plied disk warp is opposite to the one discussed in Sect j4.2.l1 
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Fig. 10. Left panel: dust continuum at 1.3 mm (color scale) and CO spirals (black arcs). Middle panel: HST image (color scale) 
and CO spirals. The 1.3 mm dust ring is marked as a white ellipse. Right panel: Comparison of 1.3 mm continuum (contours) and 
the HiCIAO polarization intensity (PI) image (color scale) at NIR. The black and red arcs mark the CO spirals identified in Fig.O 
and the HiCIAO arm (IHashimoto^ al.ll2(jllh . respectively. The SI arm is marked as a solid arc, because it is the one with apparent 
corresponding CO spirals. The coronagraphic occulting mask of (y.'3 diameter is marked as a solid- white circle. 
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Fig. 11. Left panel: plot of the normalized intensity of 1.3 mm continuum and HiCIAO PI image at NIR along the major axis of the 
disk's major axis. The central (y.'3 in HiCIAO data is occulted. Right panel: Plot of 1.3 mm intensity versus detected with the 
PdBI data of the best fit dust ring, marked as a white ellipse in Fig. [T] bottom right panel. The error bar is set to 0.23 mJy beam"^. 
The center of the ellipse is on the 1.3 mm stellar peak. The horizontal bar at the lower left corner marks the synthesized beam (Of! 45). 



above. Given the possible brightness asymmetry (see Sect j3.2.2l) , 
a better explanation is perhaps that the interpretation in terms of 
Keplerian disk is totally inappropriate. Instead, this gas may be 
tracing streamers of gas infalling (and rotating) from the ring 
onto the central object. 



4.3. CO spirals 



Our observations show that the disk structure in AB Aur is 
highly asymmetric. Although sensitive CO 3-2 and CO isotopes 
are still needed to better constrain the temperature and density on 
the spirals, our data provide the dynamics of this morphological 
component for the first time. 



4.3.1. Geometry 

Spiral-like structures have been detected both in the CO gas 
and the dust emission/scattering light at various wavelengths, al- 
though they do not come from exactly the same locations. In the 
IR, Fukagawa et aL| (120041) report four spirals (one inner arm, 
two outer arms with a branch o n one of the outer arms). More 
recently, Hashimoto et al.l (l201ll) identified eight spirals, namely 
SI, S2, S8, based on the infrared polarization-intensity im- 
age. In comparison, the CO S2 shares part of the base with the 
infrared SI (inner arm in Fukagawa's spiral) and S2 arms, which 
are in the southern part of the 1.3 mm dust ring (Fig.[TO|). The 
CO S3 arm seems to trace the inner part of infrared S3 sprial 
(Fukagawa's outer a rm). We note that the spiral detected at 5.68 
km s"^ in CO 3-2 bylLin et al. (2006) appears in between CO SI 
and CO S2 and traces the NIR SI spiral more, suggesting that it 
is at a more upper layer than the CO 2-1 spirals. 



13 



Ya-Wen Tang et aL: The circumstellar disk of AB Aurigae: evidence for envelope accretion at late stages of star formation? 



In NIR, the spirals located in the southern region appear 
shorter, suggesting an inclination effect on a geometrically flar- 
ing structure. The spirals closer to us are compressed compared 
to those located on the other side. This is also consistent with 
the NIR data, which reveal that the southern part is brighter and 
thus closer (forward scattering). To constrain the physical condi- 
tions of the spirals, high angular-resolution and high- sensitivity 
observations of ^^CO and ^^CO at multiple J states are needed. 

4.3.2. Gas kinematics 

So far, we only have the velocity information along the line of 
sight. Assumptions about the 3-D geometry are required to re- 
cover further information. Starting from the disk geometry, the 
observed velocity of gas motions with respect to the local stan- 
dard of rest (LSR), Vlsr, can be decomposed in cylinder coordi- 
nate (r, OpA,z) : 

VlSR = Vsys + [Vrot(^)cOS(^PA) + yrad(^)sin(^PA)]sin(0 

+ yz(r)cos(0, (6) 

where V^ys is the systemic velocity with respect to LSR. We 
chose 5.85 km s"^ for AB Aur in this calculation, because the 
spiral-like structures are extended. Here, Vrot is the rotation ve- 
locity (positive in counter-clockwise), while Vrad (positive for 
outward motions) is the radial velocity in the disk plane, and 
is the component perpendicular to the disk plane. The symbol 
^PA is the position angle of each location with respect to the ma- 
jor axis of the disk (-121.3° from the north), and / is the inclina- 
tion angle of the disk axis with respect to the line of sight. When 
calculating ^pa in each location of spirals, we assume /=23°, be- 
cause this best-fit model gives a Keplerian disk (see Table [3]) . 
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Fig. 12. Plot of ^PA versus ykep (^pa) as defined in Eq. [7] of 
^^CO,^^ COdisk , ^^COspirai and a purely Keplerian rotation in 
corresponding colors. The best-fit of each component is shown 
as a curve. 



For pure Keplerian rotation, Vrot equals ^JGM/r, where M 
is the enclosed mass within a radius r, and Vrad and are 0. 
Equation [6] can be rearranged as 

Jkep(^PA) = (Vlsr - Vsys) = VGMcos(^A)sin(0, (7) 
where ykep(^PA) now only depends on cos(^pa). 

Similarly, if we can correct for radial dependencies, pure ra- 
dial motions would vary as sin(^PA). Fo r examp le, for free-fall 
accretion in the disk plane, Vrad equals ^/2GM/r, Vrot and are 
0. Equation [6l can be rearranged as 

yMiOpA) = (Vlsr - V,y,) = - V2GMsin(^PA)sin(0. (8) 

Thus, a display oiy = [Vlsr - Vgys] as a function of ^pa can 
reveal the relative contributions of (Keplerian) rotation and other 
motions. We fitted 1) a pure rotation y = ai cos(^pa - a2), and 2) 
a combination of rotation and radial motions y = a^, cos(^pa) + 
sin(^PA) to the observed velocities in the spirals. We excluded 
Location IE from the fit, because the two velocity components 
are not separated well here. We also exclude CO S4 in the fitting, 
because there is less excess emission on the spirals. The fitting 
results of a pure rotation are shown in Fig.[T2]and in TableO For 
illustration, a pure Keplerian rotation with M=2.4 Mq (based on 
the spectral type of AB Aur) and /=23° (based on the axis ratio 
of the dust ring) is also shown. We note that this choice of a rep- 
resentable Keplerian rotation in AB Aur is not trivial, because 
we do not find a consistent model based on our 1.3 mm con- 
tinuum map and the high-velocity ^^CO J=2^1 line (see Sect. 
[3X2I and [3^21) . Nevertheless, the chosen M of 2.4 M© and / 
of 23° can describe the general trend of the derived y values of 
the ^^CO disk component and also the ^^CO line. We also plot a 
Keplerian rotation with M=2.4 M© but with /=-20° in Fig.Oin 
order to demonstrate that the ^^CO spiral component can originate 
in corotating gas with high inclination angle with respect to the 
disk plane. See section l44l for more discussion. 

Despite the previously noted apparent devia tions from 
Keplerian motions (best-fit velocity exponent 0.42, iPietu et aP 
2005), the disk component is very close to pure Keplerian rota- 
tion, because the best-fit value ^2 of Voco and ^12^0^1^], are 3° ±4° 
and 2° ±5°, respectively. The kinematics of the spiral component 
shows more deviations from a simple cosine function. To the 
first order, it is fitted with a similar ^2 - but with negative ve- 
locity. In other words, the kinematics of the spiral is apparently 
dominated by rotation in the opposite direction to the disk spin, 
with little contribution from radial motions. This is confirmed 
by the fit of rotation -h radial motions. The best fit gives a^, of 
-1.09 ± 0.10, which is ^ and ^4 of 0.41+0.22, which is only 
marginally diff'erent from 0. That ^4 is 2.5 times smaller than a^, 
suggests that the motions within the spirals are mainly rotation 
with slightly radial motions. 

Table 5. Best fit of ykep(^PA) 





''CO 


^^COdisk 


i2rn 

spiral 


ai 


1.36 + 0.04 


1.41 +0.05 


-1.17 + 0.12 


a2 


3 + 4 


2 + 5 


-21 + 10 



In this decomposition, motions perpendicular to the disk 
plane would not depend on ^a- The symmetric motions perpen- 
dicular to the disk midplane cannot be seen: these would only 
add up velocity dispersion after considering the limited angular 
resolution, but not change the projected mean velocity. 



4.4. Overall picture of the AB Aur system 

Gas and dust clearly indicate the existence of an inner cavity, 
of radius ^100 AU. As AB Aur is a strong UV emitter, a pos- 
sible explanation of such a cavity could be photoevaporation. 
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Fig. 13. Illustration of the origin of the spiral emissions. The 
flared disk is plotted in gray segments and arcs. The solid fil- 
aments are the detected CO spirals. The hatched filaments with 
relatively larger inclination angles are not detected. The thick 
lines mark the large-scale, flattened, and tilted infalling enve- 
lope. 



Following [Alexander et al.l (l2QQ6l) . a transitional "hole" can ap- 
pear during the evaporation of the disk. The disk mass in AB 
Aur is ~ 0.01 Mq, which is a factor of ten more massive than in 
their model. In addition, the inner disk is not completely evapo- 
rated since it has been detected at NIR. The high accretion rate 
of AB Aur (~10"^ M© yr"^) furthermore suggests that the in- 
ner disk is still replenished. This inner disk can provide substan- 
tial attenuation for ionizing radiation from the central star, and 
photoevaporation cannot be the dominant process in clearing the 
disk. Besides, the asymmetry of the dust ring has no simple ex- 
planation in this scenario. The AB Aur disk probably cannot be 
explained by such a model. 

An alternative formation mechanism of the gap/cavity is that 
the disk is in tidal interaction with an embedded companion or 
(proto-)planet. Accounting for the inner disk, the width of the 
continuum gap is ~ 90 AU. Resonant motions with a compan- 
ion in the gap can produc e high-contrast inhomogeneities (e.g. 
iKuchner & Holma nl l2003h . which would naturally explain the 
observed azimuthal brightness variations in the ring. However, 
the companion mass should be high enough so that the gap can 
re main mostly clear of gas or dust, as in the simula tion work 
of iKlev et al.1 (l200lh . Following Imeuchi et al.1 (Il996b . the half- 
width of a gap, w, opened by a low-mass companion in the 
case of a circular orbit can be estimated from: w = 1.3 r A^^^, 
where r is the semi-major axis of the companion orbit and 
A the strength ratio of tidal to viscous efl'ects given by A = 
(Mp/M^)^(3a(h/r)^)~^ . With standard values for the ratio of the 
scale height h to the orbit radius r, h/r ~ 0.1 and the viscos- 
ity parameter, a, ~ 0.01 and assuming a stellar mass, of 2.4 
Mo, a gap of half- width w of 45 AU located at r ^ 45 AU 
can be opened by a b ody of mass Mp ^ O.OSM©. However, 
[Hashimoto et al.l fcoill) rule out companions more massive than 
about five Jupiter masses from IR polarization intensity images. 
Based on an extreme assumption, 100% of polarized emission 
for the planet, this may lead to a low value of the upper limit. 
Nevertheless, having one or multiple lower mass bodies in ec- 



centric orbit or a lower viscosity (since the minimum mass scales 
as ^/a) may partially alleviate these discrepancies. 

Although the NIR spirals appear more consistent with a sim- 
ple flared disk model, we cannot explain the CO spirals with 
such a model. At NIR wavelengths, opacity prevents us from 
seeing throughout the disk, and we only see the front surface. In 
the northern part (the farside), the front surface is projected at 
a lower inclination angle, while in the southern part (the near- 
side), it is projected at a higher inclination angle. This results in 
a compression of the NIR spirals toward the south (see Fig. [13] 
for an illustration). However, the CO opacity at the spiral ve- 
locities is moderate, so we could in principle see throughout the 
disk. An explanation of the CO velocities in terms of difl'erent 
inclinations of the front and back surfaces of a flared disk would 
only be possible if we saw two anomalous velocity components 
symmetrically placed around the normal disk velocity. The CO 
spirals must originate in a structure that is more complex than a 
flared disk. 

Are the CO spirals body-driven? It has been shown numer- 
ically that a companion in a disk is expected to drive a single- 
armed spiral wake (e.g. Ogilvie & Lubow 2002). The multiple 
arcs/spiral arms observed in the NIR and mm appear difficult to 
reconcile with a single companion. Multiple bodies spread be- 
tween distances ranging from 30 to 70 AU would offer a better 
alternative. However, the pitch angle of such spirals is typically 
very small. The biggest difficulty with the body-driven spiral pat- 
tern is the apparent retrograde motion of the gas. This cannot 
be obtained with purely in-plane motions. The easiest explana- 
tion for the change of the sign in velocity is an inclination an- 
gle difference. This could relate to the indications of disk warp. 
In the multiple-body hypothesis, resonant interaction between 
these bodies may pump up orbital inclination (and eccentric- 
ity) changes, and perhaps more naturally lead to highly (> 35°) 
inclined spiral patterns. However, simulations of multiplanet -h 
disk systems have shown that inclination and eccentricity are in 
general damped quite efficiently in such cases (see review by 
iKley & NelsonI 120121 and references therein). Finally, we note 
that density contrast in companion-driven density waves is gen- 
erally significant (10% to 30 % or so), while the features we have 
detected here in CO are much fainter. 

Counter-rotating gas and spiral pattern may also trace a re- 
cent ffy-by. Such an event would na turally induce a significant 
warp in the disk (iNixon & Pringle 2010). Traces of the event are, 
however, exp ected t o decay on a timescale ^damp - ^out(2/7r/(2) 
(iLubow et al] l2002h . where Pout is the orbital period at the 
edge of the outer disk. The ^damp is on the order 10^ yr for 
a = lO'^. lPietu et al.l (l2005 l) suggest that two young stars might 
have had an encounter with AB Aur on such a timescale, JH433 
( Jones & Herbig"l979') at any time older than 35 000 years ago, 
and RW Aur some 500 000 years ago. In this respect, it is inter- 
esting to note that RW Aur is a very peculiar T Tau ri star binary 
(1.5 M q for the primary, 0.8 M© for the secondary JWoitas et al.l 
I2OOI ). with a high accretion rate, which exhibits clear evidence 
of tidal in teractions between the two stars, suggesting an eccen- 
tric orbit (ICabrit et al.ll2006h . However, the proper motion mea- 
surement accuracy of the suggested ffy-by stars precludes accu- 
rate conclusions. 

An additional important piece of information of the kine- 
matics provided by our observations is the mis-alignment be- 
tween the "large-scale" (> 2000 AU) velocity gradient in the 
cloud and the AB Aur disk axis. If interpreted by rotation, it in- 
dicates that the speciffc angular momentum of the surrounding 
envelope projects nearly perpendicular to the AB Aur disk axis. 
Material from this envelope should thus accrete with very dif- 
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ferent kinematics from that of the disk, perhaps explaining the 
apparent counter-rotating spirals. 

The large-scale velocity gradient might also be interpreted 
by infall motions. From the scattered light images, the southern 
part is closest to us, and from Fig.[3l this gas is indeed red-shifted 
as expected for infall. The CO spirals could be traces of en- 
hanced accretion along filaments. In the simplest model, which 
assumes infall from an envelope with rigid-body rotation, accret- 
ing materi al is expected to foll ow quasi-parabolic, nonintersect- 
ing orbits (ICassen & Moosman 1981). In this case, our "spirals" 
may just trace some parts of such parabolic orbits, and their ap- 
parent geometry remains consistent with such an interpretation. 
In the orbit plane, the pitch angle of the "spirals" is given by the 
ratio of radial to rotational velocity 

tan(^p,t) = S Vrad/Vrot, (9) 

where ^ is a sign that depends on conventions for velocities. This 
intrinsic pitch angle can be recovered from the "spiral" geom- 
etry by assuming some inclination angle / for the orbit plane. 
For |/| < 30°, the correction for inclination will remain small. 
The shape of our "spirals" is consistent with apparent counter- 
clockwise rotation (i.e. direct rotation in the left-handed RA,Dec 
system) and inward motions. For / ^ -20°, i.e. an orbit plane in- 
clined by about 45° compared to the main quasi- Keplerian disk 
of AB Aur at / = -i-23°, this direction of rotation is consistent 
with the overall rotation of the system. The infall motions then 
contribute to red- shifted gas along the line of sight in the north- 
ern regions (CO Spiral SI and S2) and blue-shifted gas in the 
southern part (CO S3). For S4, which lies close to the disk's ma- 
jor axis, the projection of the infall motions would essentially be 
zero. This behavior is entirely consistent with the observed ve- 
locity shifts in Fig.fSj In addition, the jkep values of the ^^CO spiral 
component follow the trend of a pure Keplerian rotation with M 
of 2.4 M© and / of -20° (Fig.[T2l), suggesting that the spiral gas 
develops from a higher inclination with respect to the disk plane. 
The analysis combining both rotation and radial motion suggests 
that the gas motion within the spirals are mainly rotation with 
slightly radial motion. A schematic is given in Fig.[T3l 

We do not see evidence of any gas coming from the opposite 
side of the main disk (dashed clumps in Fig. [13]). This may be 
the result of the combinations of two eff'ects. First, the "spirals" 
are visible because they trace density fluctuations. Such fluctu- 
ations are inherently random, so we just observe one snapshot 
of events. Second, similar features located on the opposite side 
would be much more inclined along the line of sight, around 
70°. Thus, they would only project much closer to the center of 
mass and at relatively high velocities. As infall gives a velocity 
gradient along the project minor axis, while rotation leads to a 
velocity gradient along the major axis, combination of this gas 
with the main rotating disk would lead to a distortion of the iso- 
velocity pattern, which could be difficult to identify if the amount 
of infalling gas is small as in the identified spirals (clumps in 
solid lines in Fig.[T3l). 

The accretion rate through the spiral can be estimated by di- 
viding the spiral mass by the free-fall timescale at 1000 AU; this 
gives about 10"^^ to 10"^ M© yr"^. This must be understood as 
a lower limit for the total envelope accretion rate, since gas with 
no specific overdensity is not accounted for in this process. 

Infall mo tions are conimon in Class pr otostellar envelope s, 
(e.g. L1527 lOhashi et al.1 11997'), (L1544 TafaUa et al.1 [1998), 
as suggested by the predominance of blue- shifted asymme- 
tries i n the line profiles obs erved at moderate angular resolu- 
tions (iMardones et al.l IT997b . Nonaxisymmetric accretion has 



also been invoked by lTobin et al.l (l2012h to explain the kinematic 
structure of Class protostellar envelopes seen at high angular 
resolution (Tobin et al. 201 1). The AB Aur case is special, how- 
ever, because the estimated age of the star is 1 to 4 Myr (see 
discussion in.Pietu et aD [2Q05i) . As the free-fall time scale at 10 
000 AU is on the order of 10^ years, and scales at most as r^^^ 
(when the envelope mass is negligible), material still accreting at 
this stage should come from initial distances around 50 000 AU, 
unless some mechanism has slowed down the infalling material. 
The change in magnetic sup port due to ambipol ar diffusion is a 
possible candidate (see e.g. iMouschoviasll 1 996l and references 
therein). 

AB Aur is for the time being the oldest object in which infall 
motions may have been detected. However, the apparent pecu- 
liarity of AB Aur in this respect may be a selection effect. AB 
Aur is a relatively massive star, and most likely the most massive 
of all stars with disks already studied in the Taurus clouds. Most 
studies of disks around young stars have instead been performed 
on stars isolated enough from their molecular cloud to reduce 
confusion and focus on the disk properties. Accordingly, very 
little is known about the history of accretion from the envelopes 
at large ages. 

5. Summary 

We present in this paper new mm continuum and CO observa- 
tions with an angular resolution as high as 0^.^37 (50 AU). The 
key results are summarized below. 

1. Unresolved 1.3 mm emission was detected toward the star. 
The spectral index between 3.6 cm and 1.3 mm is consistent 
with optically thick, free-free emission in a jet. However, 
the 1.3 mm emission may also trace an inner dust disk, of 
radius smaller than 50 AU, and (minimal) gas-hdust mass, 
Mgas, -4x10-5 Mo. 

2. A dust gap with a width of ~ 90 AU was then clearly de- 
tected. It can be formed with a companion mass 0.03 M© at 
a radius 45 AU. Since the outer disk is highly asymmetric 
and more massive than the disk mass predicted by the pho- 
toevaporation model by a factor of 10, our results are more 
in favor of the existence of a companion. 

3. At a radius of -145 AU from the star, a resolved dust ring 
was observed. We find that this ring is inhomogeneous in 
both radius (lower-middle panel of Fig. [T]) and azimuth (Fig. 
[TTT) . The 1.3 mm flux density is 80 mJy, corresponding to 
Mgas of 5x10"^ M©. 

4. Emission of the CO J=2^1 line was also observed within 
the dust ring extending inside the dust cavity down to a ra- 
dius < 20 AU. The axis of this CO inner disk is at a positional 
angle of -38° (Table [3]). This CO gas is likely rotating ,but 
at non-Keplerian velocities with an inclination angle of 29°, 
which is intermediate between the large-scale inclination of 
42° and the inclination angle of < 20° derived from NIR in- 
terferometric measurements. This suggests that the structure 
is warped. The origin of this inner CO disk is probably mul- 
tiple. This may be gas located in the inner part of the ring 
contaminated by possible streamers flowing onto the inner 
dust disk through the dust cavity. 

5. Four CO spirals with lengths ~ 500 AU were detected based 
on the integrated intensity and the velocity dispersion maps. 
The Mgas is 10"^ < Mgas < 10"^ M© and the kinetic energy, 
£^kin, c ontained in the spirals is 10^^ < £^kin < 10^^ erg (Sect. 
13.3.31) . In comparison to the disk mass of 0.01 M© and the 
possible companion, Mgas and Ey±^ in the spirals are very 
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small. A small perturbation of the companion with a mass 
of 0.03 M© at a radius of 45 AU can release enough energy. 
However, the multiple spirals are difficult to explain with a 
single perturbing body. Two of the CO spirals (CO S2 and 
CO S3) have their infrared counterparts (SI and S3). The 
base of the CO S2 is seen to be part of the dust ring (Fig. 

6. The kinematics of the excess ^^CO J=2^1 emission on the 
spirals were also analyzed. We found that this small amount 
of gas is apparently counter-rotating with respect to the rota- 
tion of the ring/disk system. Such behavior is difficult to ex- 
plain through planet formation scenario. A hypothetic fly-by 
(perhaps the nearby binary RW Aur) may lead to this unusual 
situation. 

7. Using the large-scale velocity gradient measured with the 
30-m, we suggest that the most likely explanation for the 
spirals and their peculiar kinematics is inhomogeneous gas 
accreting well above and/or below the midplane of the main 
disk. This explains the apparent counter-rotation by a signif- 
icant change (~ 45°) in inclination angle. It would indicate 
that weak accretion could still occur 1 to 4 Myr after the ini- 
tial star formation. 
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